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ABSTRACT 

We argue that the energy injection of pulsations may be of greater importance to the mass-loss rate of AGB stars than metallicity, and 
that the mass-loss trend with metallicity is not as simple as sometimes assumed. Using our detailed radiation hydrodynamical models 
that include dust formation, we illustrate the effects of pulsation energy on wind properties. We find that the mass-loss rate scales 
with the kinetic energy input by pulsations as long as a dust-saturated wind does not occur, and all other stellar parameters are kept 
constant. This includes the absolute abundance of condensible carbon (not bound in CO), which is more relevant than keeping the 
C/O-ratio constant when comparing stars of different metallicity. The pressure and temperature gradients in the atmospheres of stars, 
become steeper and flatter, respectively, when the metallicity is reduced, while the radius where the atmosphere becomes opaque is 
typically associated with a higher gas pressure. This effect can be compensated for by adjusting the velocity amplitude of the variable 
inner boundary (piston), which is used to simulate the effects of pulsation, to obtain models with comparable kinetic-energy input. 
Hence, it is more relevant to compare models with similar energy-injections than of similar velocity amplitude. Since there is no 
evidence for weaker pulsations in low-metallicity AGB stars, we conclude that it is unlikely that low-metallicity C-stars have lower 
mass-loss rates, than their more metal-rich counterparts with similar stellar parameters, as long as they have a comparable amount of 
condensible carbon. 
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1. Introduction 

According to the most commonly accepted scenario, the mass 
loss of AGB stars is driven by radiation pressure on dust-grains 
in the outer stellar atmospheres. This pressure causes the accel- 
eration of the outer layers, which develop, as the escape-velocity 
is reached, into a stellar wind. This mechanism is enhanced by 
the pulsations of these stars and the resulting shock waves which 
create conditions for efficient dust condensation. The wind mo- 
mentum is, however, generated by radiative acceleration, since 
pulsations alone cannot inject sufficient energy to sustain an out- 
flow (Wood 1979). 

The metallicity of an AGB star can affect its wind prop- 
erties since it influences the atmospheric structure and possi- 
bly the availability of condensible material for dust formation. 
Whether low metallicity should be associated with smaller mass 
loss is however unclear. Bowen & Willson (1991) argued that 
low-metallicity AGB stars should experience a so-called super- 
wind phase, but proposed that this should occur at a higher lu- 
minosity than for a metal-rich case. Their argument was that at 
lower metallicities the molecular opacities are smaller, which in 
turn can reduce the extension of the atmosphere. Since the on- 
set of dust formation depends on whether the atmosphere has 
reached a critical size, it appears correct to expect that weaker 
winds occur at lower metallicities. On the other hand, obser- 
vational studies of metal-poor carbon stars in the Magellanic 
Clouds and the Milky Way indicate similar mass-loss rates 
jvan Loon et al. 2000l:lGroenewegen et al. 2007h . 

We report results of an on going project to explore the mass 
loss properties of C-rich AGB stars as a function of stellar pa- 
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rameters. We question the common idea that AGB stars at lower 
metallicity should have weaker winds, and argue that the in- 
jection of kinetic energy by pulsations plays an important role 
in wind formation and the resultant mass-loss rate. We claim 
that commonly-made assumptions about boundary conditions in 
mass-loss models of AGB stars make it difficult to compare re- 
sults in the literature. The simulation of stellar pulsations, using 
a time-dependent inner boundary condition, requires, in partic- 
ular, a careful choice of parameterisation; for example, results 
derived using different rates of kinetic-energy injection, should 
not be compared. 

Theoretical studies of dust-driven AGB m ass loss at sub- 
solar (LMC) metallicity (' Wachter et al. 20081) . suggested that 
these stars should have slower winds, but that the mass-loss rate 
should be mostly unaffected by the decrease in metallicity. In 
these models, a constant C/O = 1.8 was assumed, which was 
the same value as in the solar-metallicity models of Wachter 
et al. (2002). Due to this assumption, the LMC-models have a 
significantly lower abundance of free carbon, i.e., carbon that 
is not bound to CO and may therefore form dust particles. As 
a consequence, less dust forms. A similar effect was observed 
in the results by Helling et al. (2000), who found that models 
of LMC metallicity have, in genral, slower winds and a lower 
mass-loss rate. Given the assumptions, these trends are unsur- 
prising, since reducing the amount of free carbon, and conse- 
quently the dust/gas-ratio, makes momentum transfer from radi- 
ation to dust and gas less efficient, leading to a slower wind. We 
note that molecular bands in metal-poor C-stars are as strong as, 
or stronger than the correspondi ng stars at solar metallic ity. This 
was first observed in the optical (IWesterlund et al. 199 lb and re- 
cently confirmed in the infrared (see, e.g., Lagadec et al. 20071 
and refereces therein). The wind models presented by Wachter 
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et al. (2008) and Helling et al. (2000) are also based on grey 
radiative trans fer, which appears to affect the mass-loss rates 
(author?) (cf. iHofner et al. 2OO3I) . Moreover, the velocity am- 
plitudes and not the rates of energy-injection by pulsations were 
used to quantify the effects of pulsations on mass loss (see, also, 
[Winters et al. 20001) . 

2. Numerical Modelling 

The models for the atmospheres and winds are computed us- 
ing frequency dependent radiative transfer for the gas and dust, 
including detailed micro-ph ysics of th e dust gr ains and their for - 
mation (described in A ndersen et al. 2003: Hofner et al. 2003h . 
in combination with time-dependent hydrodynamics. Here we 
describe the initial and boundary conditions, in addition, to fur- 
ther assumptions. 

2.1. Pulsations and the Inner Boundary Condition 

The kappa-mechanism (see, e.g.. l01ivier & Wood 20051) . which 
is believed to be the cause of long-period pulsations, operates in 
the convective envelope. In our models, we simulate the effect of 
this mechanism on the atmosphere, using a time-dependent inner 
boundary condition, i.e., a so-called piston boundary of the form 

Rinit) = RiniO) + £amp Sin J^fj , (1) 

where Awp is the velocity amplitude and P is the pulsation 
period. We emp loy an empirical period-luminosity relation 
dFeast et al. 1989h . which maintains a direct relationship be- 
tween the period and luminosity, for all models. Since the bottom 
layers of the atmosphere must be optically thick for the model 
atmosphere to be realistic, one should place the piston boundary 
in that optically thick region of the atmosphere. This leads to in- 
ner boundary conditions for the metal-poor (SMC-like) models 
where the gas pressure is significantly higher at R - Ri„ (see 
Fig.[T]), while the temperature changes much less. Thus, the ki- 
netic energy-injection by the pulsations/piston is typically larger 
in these models for a given Au^. 

2.2. Free Carbon and Dust Formation 

The raw-material for dust grains is the condensible carbon in the 
atmosphere, i.e., the carbon that is not bound in CO molecules 
and therefore, in principle, available for dust formation. Since 
the absolute oxygen abundance so may vary depending on the 
initial metallicity of the star and the amount of oxygen dredge- 
up, sc ~ sc- So is a more relevant quantity than C/O = sq/so- 
Low-metallicity C-stars do not necessarily have carbon-poor 
atmospheres. Since oxygen (along with large amounts of car- 
bon) isjiwstlikely dredged-up in each thermal-pulse cycle (see, 
e.g.. lHerwig 20051 and references therein) the C/O-ratios do not 
have to be large, even if ec is relatively large. For C-stars with a 
solar (or higher) relative abundance of oxygen, the oxygen abun- 
dance is hardly affected by oxygen dredge-up, but for metal-poor 
C-stars, where the inital abundance is small, it certainly is. We 
emphasise that observations by Wahlin et al. (2006) indicate that 
C/O-ratios in metal-poor C-stars may be almost an order of mag- 
nitude higher than at solar metallicity; so even without oxygen 
dredge-up, ec may therefore be quite large. After a number of 
thermal pulses, there may be as much free carbon in a metal- 
poor C-star atmosphere as for a solar-metallicity. 



Table 1. CNO abundances and metallicity of the models con- 
sidered in this letter. 





log(£c) 


log(£c) log(eN) 


log(eo) [Fe/H] 


AGS05 


8.50/8.80 


8.89/9.04 


7.78 


8.66 0.0 


Sub-solar 


8.50/8.80 


8.81/8.86 


6.78 


7.96 -1.0 


Table 2. 


Parameters 


used for the models. 


All models have 


M^, = IM0, A-models have log(e) 


= 8.80, and B-models 


log(e) = 8 


.50. 










Model 


log(L/Lo) 


7-eff (K) 


log(g) 




Aifp (km s ' ) 


Al 


3.85 


2800 


-0.65 


390 


3.0 


A2 


3.85 


2800 


-0.65 


390 


4.0 


A3 


3.85 


2800 


-0.65 


390 


5.0 


Bl 


3.70 


2600 


-0.64 


294 


4.0 


B2 


3.70 


2600 


-0.64 


294 


5.0 


B3 


3.70 


2600 


-0.64 


294 


6.0 



We assume that s is almost unrelated to the over all ini- 
tial metallicity of the star. We use the same s for subsolar and 
solar-metallicity models, i.e., we compare the effects of different 
metallicities only. Since e, and not the C/O-ratio, remains fixed, 
it is not obvious that a lower initial metallicity leads to lower 
mass-loss rates, when comparing mass-loss models at different 
metallicities. 



3. Results and Discussion 

3. 1 . Numerical Results 

We computed a set of models, where we considered two sets 
of stellar parameters (labelled A and B), two chemical composi- 
tions (see Table[T]and Table|2]i, and various positions of the inner 
(piston) boundary. The results are the following (see also Fig.|2] 
and Fig. [3). For a fixed value of ec and a given set of stellar 
parameters, we found that: 

- The location of the inner boundary can affect significantly 
the mass-loss rate if Aup is kept fixed. 

- Models with a metal-poor composition (apart from carbon) 
have similar, or higher, mass-loss rates compared to models 
of solar composition (defined as that of Asplund et al. 2005) 
if the kinetic-energy injections by pulsations are similar 

- Models with a solar composition, in general, show slower 
winds, lower mass-loss rates, and lower degrees of dust con- 
densation than models with metal-poor composition for fixed 
Awp if the inner boundary is located at smiliar optical depth, 
due to the difference in density at that point. 

We emphasise that we compare models with the same absolute 
abundance of free carbon, which excludes the possibility that 
these effects are due to the accessibility of raw material for dust 
formation (see Sect. 12.21 for further details). However, the dif- 
ferent absolute abundances of elements other than carbon affect 
molecular opacities and, consequently, the structure of the atmo- 
sphere. Compared to these abundance differences, the diflFerence 
in the mass-loss properties of both of the models considered (A 
and B), appears to be more strongly correlated with the input 
of kinetic energy at the inner boundary. In our models, we use 
a sinusoidal piston boundary condition with a velocity ampli- 
tude AMp, which means that the energy injected into a mass shell 
of mass dm at the inner boundary is dEp = ^Au^dm. From an 
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Fig. 1. The Rosseland mean opacity at the in- 
ner boundary of the hydrostatic initial atmo- 
spheric structures as function of the geometric 
location of the boundary and the gas density 
at the boundary (at / = 0) for solar (Asplund 
et al. 2005) and a typical metal-poor compo- 
sition ([Fe/H] = -1.0). The full line in the 
left panel shows an opacity-law of the form 
where tq is the opacity at 



TRo, 



Toe 



r/Ro + l 



a some radial point Rq- The horizontal dashed 
lines marks tross = 1. Points below this line 
are usually associated with too shallow initial 
structures. 



observational point of view, the mass loss rate appears to be pro- 
portional to the infe rred pulsational energy injection rate of AGB 
stars (see Ivan Loon 2008. and references therein). We therefore 
define a quantity 

q = Aulpinit = 0), (2) 

which is proportional to the kinetic-energy injection, and com- 
pare with the wind properties. We find that, in general, 

log(M) cc \og(q), (3) 

while the wind velocity Mqui shows a more complicated depen- 
dence on q (seeFig.|2]and Fig. [3]). The latter is caused by the levi- 
tation of the outer atmospheric layers and the increase in density, 
which affects the dust-formation efficiency. 

In Fig. [2] and Fig.|3] we note that the B -models with metal- 
poor composition, show typically faster winds than solar metal- 
licity B-models, and that metal-poor A-models show mass-loss 
rates higher than the corresponding models with solar compo- 
sition (see Fig. |2]l. We note also that the mass-loss rates and 
mean degrees of dust condensation (/c-values), obtained for the 
A-models with metal-poor composition, in principle, show no 
trends with q; in these cases, this is because the region of wind 
formation is saturated with dust, i.e., /c reaches its upper limit 
(typically ~ 0.4 - 0.5), and the wind is maximally dust-driven. 
Furthermore, the metal -poor A-models with AMp $ 2kms"' can- 
not sustain a wind due to inefficient dust formation. The same is 
true for B-models with AMp S 3 km s"' . This is because the tem- 
perature gradient is slightly flatter in the AGS05 models, which 
affects the conditions for dust formation negatively and therefore 
requires a higher kinetic-energy input to trigger dust formation. 

3.2. The M-q Connection 

From the above, it is obvious that there is a connection between 
the average mass-loss rate and the <7-parameter. In principle we 
find three different cases: M oc q^^^, M cc q, and M - const., 
which can be derived from the equation of motion for the atmo- 
sphere and wind; a detailed discussion is beyond the scope of 
this letter and will be published in a forthcoming paper. Before 
we discuss the consequences (see Sect. 13.31 ), it is of interest to 
discuss how these dependences on q arise from a physical point 
of view. 

The first two cases reflect the effects of the energy-input by 
pulsations, which affects the levitation of the outer atmospheric 
layers and mean degree of dust condensation. Hence, the mass- 
loss rate increases as more kinetic energy is pumped into the at- 



mosphere. The third trend, or the absence of a trend, is due to sat- 
uration of the dust-formation mechanism, i.e., there is an upper 
limit for the dust/gas-ratio. When this limit is reached, the dust 
opacity does not change when more kinetic energy is pumped 
into the atmosphere and the mass-loss rate remains constant. The 
transitions between these cases depend on the period/luminosity 
of the model stars, and their atmospheric structures. 

3.3. Stronger Winds at Lower Metaiiicity? 

As is evident in Fig. |2] and Fig. [3] a larger ^-value corresponds, 
in general, to a larger mass loss and the models with a metal- 
poor composition have typically faster winds and higher mass- 
loss rates for given stellar parameters. They also require larger 
^-values to sustain a wind. In the case of the A-models, the 
trends (mass-loss rates as a function of q) appear to be differ- 
ent for the two metallicities considered, although the number 
of models is too small to establish the existence of two differ- 
ent trends statistically. However, in all cases, our modelling re- 
sults suggest that the winds may become stronger at sub-solar 
metallicities, or at least not weaker than at solar metaiiicity, pro- 
vided that the abundance of free carbon is the same even if the 
over all metaiiicity is different. This result may explain obser- 
vations of mass-los s rates of C-stars in the LMC, and other lo- 
cal dwarf galaxies (van Loon et al. 20001: Ivan Loon et al. 20031: 
iMatsuura et al. 2007 1). in which the mass-loss rates appear to be 
as high, or even higher, than for C-stars in the solar neighbour- 
hood. Similar results were obtained by Jackson et al. (2007a; 
2007b) for C-stars in other dwarf galaxies. 



4. Summary and Conclusions 

The results presented above may be summarised as follows: 

The mass loss of C-rich AGB-stars, for a given set of stellar 
parameters and abundance of free carbon, is a function of the 
energy-injection by pulsations as long as the dust formation has 
not reached the saturation limit. When this limit is reached, the 
mass-loss rate no longer increases with the energy-injection. 

The range of realistic, or even possible, pulsation energies is 
limited, and connected to the stellar parameters. The treatment of 
the inner boundary and pulsations in the case of a non-saturated 
wind is crucial. Hence, the importance of constraints on the 
energy injection by pulsations cannot be emphasised enough. 
A quantitatively-coiTect theory, or reliable observational con- 



4 



L. Mattsson et al.: Models of mass loss from C-rich AGB stars at different metallicities 




-0.5 0,0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 

log(q) log(q) log(q) 



Fig. 2. Comparision between mass-loss rate, mean degree of dust condensation (/c), and wind velocity for A-models, where log(e) = 
8.8 (see also Table |2|l, as functions of the energy-injection parameter q for solar (Asplund et al. 2005) and a typical metal-poor 
composition ([Fe/H] - -1.0). The solid and dashed lines in the left panel indicate (M) oc q and (M) = const., respectively. 
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Fig. 3. Comparision between mass-loss rate, mean degree of dust condensation (/c), and wind velocity for B-models, where log(e) = 
8.5 (see also Table |2|, as functions of the energy-injection parameter q for solar (Asplund et al. 2005) and a typical metal-poor 
composition ([Fe/H] — -1.0). The solid and dashed lines in the left panel indicate (M) oc q'^l^ and (M) oc q , respectively. 



straints, of long-period stellar pulsation are required to obtain 
a complete quantitative understanding of mass loss on the AGB. 
Since the results here are based on only a rather small grid of 
models, further study of boundary-condition effects is needed 
and will be presented in a forthcoming paper. 

However, one firm conclusion can be drawn from present re- 
sults: low metallicity should not necesserily be associated with 
a lower mass-loss rate. More precisely, the development of a 
strong super-wind phase is expected for C-stars at low (sub- 
solar) metallicity. The duration of the AGB -phase may, in partic- 
ular, be aff'ected and the effects on nucleosynthesis results (stellar 
yields) may be profound. 
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